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Abstract 

Ragone  plots  relating  power-density  to  achievable  energy-density  have  been  used  for  many  years  as  an  empirical  basis  for  comparative 
performance  evaluation  of  various  battery  systems.  However,  their  theoretical  basis,  and  the  relative  contributions  of  ohmic  and  activation 
(Tafel)  polarization  effects  to  the  dependence  of  energy-density  on  operating  power-density  have  been  little  studied.  Quantitative  modeling 
calculations  for  various  selected  values  of  ohmic  resistance,  and  Tafel  exchange -current  densities  and  polarization  slopes,  b,  provide  ration¬ 
alizations  for  the  shape  of  Ragone  diagrams  and  their  positions  in  the  log-log  plane  of  energy-density  versus  power-density  plots.  Separate 
evaluations  of  the  ohmic  and  Tafel  polarization  ‘components*  of  the  Ragone  relations  are  of  practical  value  in  the  evaluation  of  the  battery 
performance.  The  Tafel  value  b  is  of  special  importance.  Complications  which  arise  in  trying  to  represent  electrochemical  capacitor  behaviour 
in  terms  of  Ragone-type  plots  are  analysed;  they  are  important  for  the  design  and  operation  of  hybrid  capacitor/battery  systems  for  electric- 
vehicle  drive  trains. 
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1.  Introduction 

During  the  past  thirty  years  relations  between  specific 
power-density  (PD,  W  kg  ')  and  energy-density  (ED,  Wh 
kg1)  have  been  empirically  plotted,  often  on  logarithmic 
scales.  Such  plots  originated  in  an  historically  significant 
paper  by  Ragone  [  1  ]  at  the  May  1968  meeting  of  the  Society 
of  Automotive  Engineers  at  Detroit,  MI.  where  the  first  exam¬ 
ples  were  shown.  We  specially  mention  this  since  Ragone 
plots  are  often  found  in  battery  literature  but  the  original 
attribution  is  rarely  quoted  nor  are  interpretations  of  the  plots 
usually  given.  In  Ragone’s  paper,  the  behaviour  of  battery 
systems  for  electrically  powered  vehicles  was  already  con¬ 
sidered,  and  is  now  a  matter  of  major  importance. 

Ragone  plots  for  batteries  are  now  complemented  by  those 
for  electrochemical  capacitors  [2,3  J  especially  as  the  latter 
are  perceived  as  energy-storage  systems  capable  of  high 
power  delivery  and  high  power-level  recharging.  Of  special 
importance  is  the  hybrid  combination  of  a  high-power  elec¬ 
trochemical  capacitor  with  a  high  ED  and  high  charge-capac¬ 
ity  battery  (4],  Hence  the  PD  and  the  ED  components  of  the 
performance  of  hybrid  systems  are  important  to  have  quan- 
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titatively  evaluated  in  order  for  suitable  matching  to  be 
achieved  and  understood  on  a  fundamental  basis. 

Ragone  plots  usually  have  a  hooked  shape  corresponding 
to  a  fall-off  of  ED  as  PD  (i.e.  power  drain)  is  increased  in 
discharge  of  the  battery.  A  theoretical  basis  for  this  shape  was 
not  provided  in  Ragone’s  paper  [  1  ]  but  the  general  origin  of 
it  was  recognized  qualitatively  [I]  in  terms  of  increasing  cell 
polarization  (hence  diminution  of  cell  voltage  which  deter¬ 
mines  ED)  as  higher  power  demands  are  made  on  a  battery 
or  capacitor  at  higher  rates  of  discharge. 

A  comparison  between  Ragone  plots  for  advanced  electro¬ 
chemical  capacitors  and  Ni  /Cd  and  lead/ acid  batteries  ( spec- 
ificaiions  not  given)  was  made  by  Idaho  National 
Engineering  Laboratory  ( INEL,  9-95 ) ,  by  Murphy  et  al.  [  3  ] , 
based  on  1995-data,  see  Fig.  1 .  Quite  high  ED  in  the  7-8  Wh 
kg’ 1  range  for  3  V  non-aqueous  solvent  carbon  double-iayer- 
lype  capacitors  and  about  1-2  Wh  kg-1  for  1  V  aqueous 
capacitors  are  attainable.  Pinnacle  Research  Institute  capac¬ 
itors  ( probably  of  the  Ru02  mixed  oxide  type)  provide  about 
0.4  Wh  kg  1  at  about  400  W  kg1.  The  3  V  organic  electro¬ 
lyte  solution  capacitors  provide  PD’S  between  200  W  and  3 
kW  kg-1.  The  US  Department  of  Energy  projected  goal  is 
about  1.7  kW  kg’ Such  PD  levels  are  observed  notwith¬ 
standing  the  intrinsic  problem  with  high-area  double-layer 
capacitors  of  distributed  internal  pore  resistance  [5,6]  within 
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Fig.  I .  Ragone  plot  for  uhracapacitors  and  batteries  adapted  from  Ref.  [1,3]. 


the  electrode  matrices  and  a  consequent  range  of  resistance- 
capacitance  (RC)  time  constants  leading  to  non-uniform 
power  capability. 


2.  Theory 

2. 1.  Formal  energy-density  and  power-density  definitions 

Energy  and  power  of  an  electrochemical  power  source  are 
defined  by 

Energy  ( E)  =  Q(  A  V) ,  for  a  battery  (la) 

or 

Energy  ( E)  —  ^  Q(  A  V)  =  ^  C(  A  V) 2  for  a  capacitor  (lb) 

and 

Power  (P)  =  /(AF)  (2) 

where  A  V  is  the  voltage  established  between  the  electrodes 
of  a  cell,  and  Q  is  the  charge  stored  at  the  potential  difference, 
A  V.  If  Q  is  expressed  per  unit  weight  then  E  is  the  ED,  or  if 
Q  is  expressed  per  unit  volume,  E  is  the  volume  ED.  In  Eq. 
(2);  /  is  the  current  drawn  on  discharge.  In  practice,  the 
withdrawable  charge  is  usually  diminished  at  elevated  power 
drains  due  to  non-uniformity  of  discharge  and/or  isolation  of 
some  of  the  active  material  at  high  /  values. 

From  Eqs.  ( 1 )  and  ( 2 )  the  ratio  of  energy  to  power  drawn 
at  a  given  A  V  for  a  battery  is  Qll,  and  for  a  capacitor  is 
Q/2I,  both  of  which  have  units  of  time. 

P  can  also  been  expressed  in  terms  of  the  current  density, 
i,  defined  by  i'  =  //A,  where  A  is  the  accessible  surface  area 
of  a  given  electrode.  In  practice,  most  battery  or  capacitor 
electrodes  are  three-dimensional  porous  structures  whose  A 
values  are  usually  not  precisely  known.  Then,  it  is  useful  to 
express  the  current  density  as  a  current  per  unit  weight  (or 


volume)  of  the  active  materials,  which  gives  immediately  the 
PD  per  kg  or  1.  We  also  use  the  symbol  i  for  this  current 
density. 

For  a  battery  system  Q  is  determined,  in  an  elementary 
way,  by  Faraday  s  laws  according  to  the  weight,  W,  of  the 
active  materia;  and  the  number  of  electrons,  z ,  per  mole 
involved  in  it"  oxidation  or  reduction  and,  thus,  its  equivalent 
weight,  M/z,  where  M  is  its  molar  weight.  Thus,  for  a  sample 
weight  W 

Q = F\Y/  ( Ml  z)~  zFW/M  (3) 

where  F  is  the  Faraday  constant,  96  487  C  mol ' 

2.2.  Dependence  of  energy-density  and  power-density  on 
current  drain  for  a  battery 

Both  ED  and  PD  of  a  battery  are  determined  by  the  avail¬ 
able  operating  voltage,  AV,  in  Eqs.  (la)  and  (2),  respec¬ 
tively.  Upon  discharge  at  i,  A  V  is  diminished  jointly  by  the 
ohmic  polarization  i/?  and  the  kinetic  or  activation  (Tafel) 
polarization  r)^  =  a  +  b  In  i.  R  is  the  ohmic  resistance  of  the 
electrolyte  solution  combined  with  that  of  the  electrode 
structures. 

The  Tafel  polarization  can  also  be  expressed  in  terms  of 
the  exchange  current  density,  i0.  passing  reversibly  at  7ja  =  0 
17] 

t(i?a)  =  »nexp[r)a/fcl  (4) 

or 

In  i(  7ja)  =  In  i0  4-  r\J b  (5) 

where  b  is  the  so-called  Tafel  slope  d7ja/d  In  i,  also  defined 
as  b  =  RT/aF  where  a  is  the  transfer  coefficient  equal  to 
0.5  ±  0.05  V  for  most  simple,  one-electron  processes  [7] . 

If  there  were  no  polarization  losses,  obviously,  from  the 
above  definition  of  power  in  Eq.  (2),  power  would  increase 
linearly  with  i.  Since,  however,  decreases  of  the  thermody¬ 
namically  available  (reversible)  potentials  at  the  cell  elec- 


W.G.  Pell.  B.E.  Cvn-:uy  /  Journal  of  Power  Sources  63(1996)  255-266 


257 


trodes  arise  as  i  increases,  and  ohmic  iR  losses  arise 
simultaneously,  the  power  must  be  expressed  as 

Power  (P)  =  i(AV)-i'[AVre,-rR-E[7)u(i)Jl 

(6) 

=  i&  Vrev  -  i2R  -  iZb  !n  i  +  iZb  In  rft 

where  A  Vrcv  is  the  reversible  cell  voltage,  b's  are  for  the 
respective  anode  and  cathode  reactions  and  i0’s  have  been 
introduced  from  Eqs.  (4)  or  (5).  Eq.  (6)  involves  the  tran¬ 
scendental  functions  b  In  i  and  i  b  In  i0,  together  with  the 
linear  terms  in  *,  and  hence  can  only  be  numerically  solved. 

Since  RT/i{)F  corresponds  to  a  faradaic  reaction  resistance 
/?F  [7],  iQ-RT/ RrF,  P  can  also  be  expressed  in  terms  of  RF 


and  with 

V=EIQ  =  A  Vrev  -  iR  -Lia  +  b  In  i) 


dlnP  E  V-iR-Zb 
d  In  £  PQ/H-iR-Zb ) 
E 

Q{-iR-Zb ) 

V  ,  A 


(12) 


-iR-Zb 

From  Eq.  ( 12),  PD  increases  with  ED  for  the  conditions 


P  =  i( A  Vtcv -iR-b  \n[RFF/RT\  -b  Ini)  ( 7 ) 

where  l  lRF  measures  the  kinetic  reversibility  of  the  reaction 
and  can  he  directly  determined  by  means  of  impedance  spec¬ 
troscopy  or  from  linear  polarization  slopes  for  t}a  evaluated 
near  the  reversible  potential  [7] . 

In  Eq.  ( 6)  /A  Vrev  is  the  ideal  power  available  with  increas¬ 
ing  /;  if  AVrev  were  uninfluenced  by  the  magnitude  of  i,  P 
would  increase  linearly  with  i.  This  is  actually  observed  for 
sufficiently  low  i  values  that  the  ohmic  andTafel  polarizations 
remain  small.  From  Eq.  (6)  it  is  clear  that  although  P  will 
initially  tend  to  increase  with  current  drain  i,  it  must  eventu¬ 
ally  decrease  due  to  the  negative  terms  in  ohmic  and  kinetic 
polarizations.  Hence  there  is  a  maximum  in  power  withdraw¬ 
able  from  a  cell  with  increasing  i  determined  by  differenti¬ 
ating  P  with  respect  to  i\  thus 

API di  =  A  VKV  —  2iR  —(i) Zb/ i  -  Zb  ln( 0  +  Zb  ln( t0)  (8) 


7  jn  C>0  for  V< iR  +  Zb  (13a) 

d  In  E 

while  it  decreases  with  increasing  ED  for 

~~~<0  for  V>iR+Zb  (I3b) 

d  In  E 

and  the  maximum  in  PD  with  respect  to  ED  occurs  when 

V=iR+Zb  (14) 

i.e.  V  for  increases  with  R  and  b  values  but  P  itself 
decreases. 

Note  that  at  sufficiently  high  power  drain  ( virtually  a  short- 
circuit),  V  will  be  reduced  to  zero  so  that  iV=  V=0  which 
arises  when 

/( A  Vrcv)  =  ?R  +  iZb  In  i  -  iZb  In  t0  (15) 


=  &VKV  —  2iR ~  Zb  —  Zb  ln(0  +  Zb  ln(*o)  (9) 


-  0  for  the  maximum  in  P 


The  PUMX  arises  at  a  current  determined  by  the  transcen¬ 
dental  equation 

2iR  +  Zb  In  i  =  A Vm- +  Zb  In  i0 -Zb  (10) 


for  which  no  explicit  solution  is  possible  so  that  the  i  at  which 
Pmiul  arises  must  be  determined  numerically. 

Note  that  P ^  is  determined  both  by  i  and  In  i;  thus,  with 
increasing  i,  it  is  seen  that  the  rate  of  decrease  of  P  is  actually 
mainly  determined  by  the  ohmic  iR  term  which  is  therefore 
practically  a  more  important  factor  than  the  kinetic  polari¬ 
zation  at  high  i  unless  R  is  relatively  negligible.  However, 
there  is  additionally  the  ‘constant’  factor,  Zb,  related  to  sums 
of  the  Tafel  slopes  for  anodic  and  cathodic  half-cell  reactions 
of  the  cell.  Also,  of  course,  better  maximum  power  is  attained 
when  the  reversible  cell  voltage,  AVrov,  is  larger,  i.e.  for  a 
bigger  open-circuit  e.m.f.  or  ‘driving  force’  of  the  cell. 

Correspondingly,  a  maximum  in  In  P  arises  in  the  ln-ln 
Ragone  plot  when  d  In  P/d  In  £=0.  Thus  with 


d  In  P_/£\  d  P 
d  In  £  \P/d  E 


(ID 


or 

A  Vrcv  =  iR  +  Zb  In  i-Zb In  i0 

(16) 

i.e. 

AVrc¥  =  i£  +  7k(r) 

(17) 

Practically,  this  arises  when  i  tends  to  be  sufficiently  large 
that  the  iR  drop  virtually  completely  opposes  the  A  Vlcv  so 
that  no  useful  power  is  any  more  available. 

The  ED  ( Eq.  ( 1 )  )  also,  of  course,  depends  on  the  dimi¬ 
nution  of  AVrev  due  to  the  same  polarization  effects  that 
influence  P;  thus 

E=Q(AVw-iR-Zb\ni+Zb\niQ)  (18) 

cf.  Eq.  (6),  where  the  term  in  brackets,  equal  to  EtQ  is 
identical  with  the  term  PI  I  of  Eq.  (6). 

It  is  the  relation  between  £  in  Eq.  ( 18)  and  P  in  Eq.  (6) 
that  is  the  basis  of  Ragone  plots  for  battery  systems. 

An  additional  factor  can  sometimes  be  the  concentration 
polarization  in  certain  cells,  e.g.  with  soluble  redox  reagents 
or  fuel  cells  with  soluble  fuel,  e.g.  CH3OH.  That  is,  at  suffi¬ 
ciently  high  current  drains,  i  ~  it,  power  becomes  limited  by 
mass  transport  of  the  reagent  or  1  imiting  by  zero  concentration 
of  the  reagent  at  the  reactive  electrode  surface. 
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Concentration  polarization,  t}c,  is  defined  by  the  equation 


Here,  i|  is  the  limiting  current  that  arises  when  the  concen¬ 
tration  of  the  reactive  species  at  the  electrode  surface  is 
reduced  to  zero  corresponding  to  maximized  concentration 
gradient  at  the  electrode  surface.  In  this  case,  the  cell  potential 
varies  as 

V=A  V^-iR-UvJ-Liyc)  (20) 

where  r)a  is  the  activation  polarization  and  7)c  is  the  concen¬ 
tration  polarization. 

2.3.  Ragone  plots  for  electrochemical  capacitor  behaviour 

Ragone  plots  have  taken  on  a  new  importance  in  the  tech¬ 
nology  of  combination  of  electrochemical  capacitors  having 
high  PD  (but  relatively  low  ED)  with  high  ED  aqueous 
electrolyte  batteries  for  electric-vehicle  drive  trains  [4].  For 
such  systems  power  level  and  energy  balancing  in  the  cou¬ 
pling  between  these  two  kinds  of  energy-storage  systems  is 
of  major  importance  [4]  and  consequently  requires  quanti¬ 
tative  knowledge  of  the  Ragone  relations  for  both  types  of 
system. 

As  we  show  below,  the  factors  that  determine  the  ED 
versus  PD  relations  for  electrochemical  capacitors  are  signif¬ 
icantly  different  from  those  for  batteries  because,  on  dis¬ 
charge,  a  continuous  decline  of  output  voltage  is  an  intrinsic 
property  of  a  charged  capacitor  (determined  by  the  relation 
C=Q/V)  while,  ideally,  a  battery  discharge  voltage  can 
remain  constant  so  long  as  the  chemical  redox  species  remain 
in  a  two-phase  electrochemical  equilibrium  prior  to  exhaus¬ 
tion  of  one  or  other  of  the  reactants,  as  is  well  exemplified  by 
the  case  of  discharge  of  Li/SOCl2  primary  batteries.  Practi¬ 
cally,  of  course,  most  secondary  batteries  do  show  to  some, 
often  significant  degree,  a  decline  of  cell  voltage  upon  dis¬ 
charge  ( e.g.  with  Li  intercalation  systems,  see  Ref.  [81),  and 
this  has  to  be  taken  into  account  in  the  evaluation  of  battery/ 
capacitor  hybrid-power  systems. 

We  now  examine  how  rating  of  PD  of  a  capacitor  can  be 
treated  and  experimentally  examined. 

Formally,  there  is  an  important  difference  between  the 
principle  involved  in  rating  power  drain  from  an  ideal  battery 
and  that  applicable  to  capacitors,  defining  an  id«*:u  battery  as 
one  for  which  the  cell  voltage  through  the  discharge  half¬ 
cycle  remains  almost  constant.  As  emphasized  earlier,  an 
essentia]  feature  of  the  electrical  behaviour  of  a  capacitor  is 
a  decline  of  voltage  upon  discharge  in  proportion  (when  its 
capacitance  is  independent  of  potential  difference  between 
its  'plates’ )  to  charge  removed  through  an  external  load,  i.e. 
V—  Q! C  or  dVldQ  =  C~ l.  Hence,  withdrawal  of  charge  at  a 
rate  i  causes  a  corresponding  rate  of  decline  in  V  so  that  a 
formal  power  rating  at  a  particular  /  depends  on  the  state  of 
discharge  of  the  capacitor. 


With  the  usual  definition  of  PD  as  i V  we  have 
P-i(QIC)  or,  differentially  P=  VC{dV/dt)  (21) 

for  a  given  constant  discharge  i‘(  =dQ/df=C  dV/dt)  and, 
for  C  invariant  with  potential,  dV/dt  =  constant  for  a  given 
fixed  r. 

In  another  way,  we  can  write  the  charge  remaining  on  the 
capacitor  at  time  t,  beyond  some  initial  condition  (i)  when 
the  charge  was  Qh  as 

Q,=Q,-jiit  (22) 

The  corresponding  potential  V,  is  QjC,  i.e. 

V,=  (S,- jidt)/C  (23) 

so  that  the  PD  at  time  f,  after  commencement  of  discharge  at 
i  A  cm  “  2,  will  be 

P,  =  i(Q,~  fidi)/C  (24) 

0 

The  initial  PD  (P,)  was  iV,  or  iQj C  so  that 

P, =iV,—ijidt/C  (25) 

o 

/ 

”«"(V,— fidf/C)  (26) 

0 

which  eventually,  of  course,  becomes  zero  when 

i 

jidt/C=Q,/C=Vi 

0 

Alternatively,  a  capacitor  device  could  be  dischargedunder 
conditions  of  constant  power,  i'V= constant.  For  such  a  situ¬ 
ation,  as  charge  was  being  withdrawn  and  V  declining,  it 
would  be  necessary  to  arrange  instrumentally  that  i  would  be 
increasing  in  order  to  maintain  the  product  iV—  constant,  ■=  k 
say.  In  fact,  i  =  kV  ~ 1  so  that  the  condition 

di/dV=  —kV  ~2  (27) 

would  be  required.  V(  =  V,)  at  any  timer  during  this  constant 
power  discharge  would  be  V,  =  QJC  in  terms  of  the  charge 

Q,  remaining  at  time  t  into  the  discharge.  Hence 

di/dV=-k/(QJC)2  (28) 

i.e.  i  must  increase  as  V  decreases,  which  is  the  basis  of  this 
case. 

These  relations  show  th-t  PD  evaluation  for  a  capacitor  is 
not  as  simple  as  for  an  ideal  battery. 

Additionally,  as  for  a  battery,  the  PD  will  decline  (at  a 
given  degree  of  charge  g,held  by  the  capacitor)  with  increas¬ 
ing  discharge  rate  i  when  there  are  ohmic  polarization  losses. 
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This  will  certainly  be  the  case  with  porous,  high-area  capac¬ 
itor  electrodes,  especially  when  less  conductive  non-aqueous 
electrolytes  are  employed  as  in  the  case  of  the  higher 
operating  voltage  double-layer  capacitors.  Note  that  for 
capacitors,  there  are  normally  no  kinetic  Tatel  polarization 
losses  involved  in  the  PD  evaluation.  However,  for  electro¬ 
chemical  capacitors  utilizing  redox  or  adsorption  iseudo- 
capacitance  [9],  additional  faradaic  (Tafel-type)  polariza¬ 
tion  losses  will  arise,  as  was  discussed  earlier  for  the  case  of 
PD  versus  ED  Ragone  plots  for  batteries. 

As  with  battery  systems,  the  rated  PD  of  capacitors  will 
depend  on  the  current  drain  at  which  the  rating  is  evaluated; 
it  will  also  depend  on  the  state-of-charge  ( SOC)  in  the  capac¬ 
itor  case  much  more  than  for  most  battery  systems  except  for 
the  case  of  Li  intercalation  types  where  quasi-  or  pseudo¬ 
capacitor  behaviour  can  arise  [8], 

A  final  difference  between  capacitor-  and  battery-power 
rating  is  the  relation  between  PD  and  ED  for  these  two  types 
of  electrical  energy  storage  system.  Eqs.  (1)  and  (2)  show 
that  there  is  an  essentially  direct  relation  between  PD  (iV) 
and  ED  (QV or  0.50V)  for  a  given  initial  cell  voltage  since, 
for  both  cases,  the  factor  by  which  the  cell  voltage  is  dimin¬ 
ished  by  polarization  effects  with  increasing  i  is,  at  least 
theoretically,  the  same.  However,  in  the  case  of  capacitors, 
the  density  of  stored  energy  scales  with  the  square  of  the 
voltage  attained  at  the  end  of  charge  while  the  PD  remains 
linear  in  that  voltage,  though  for  both  properties  the  voltage 
can  be  controllably  raised  (unlike  for  the  case  of  the  ther¬ 
modynamic  voltage  of  most  battery  cells)  up  to,  or  just  below, 
the  electrolyte  solution  decomposition  potential  limit  (i.e. 
about  1 .3  V  for  aqueous  solutions  or  about  3.5-4.0  V  for  a 
non-aqueous  aprotic  solvent  solutions) .  The  above  statement 
implies  the  interesting  conclusion  that  the  ratio  of  ED  to  PD 
for  a  given  capacitor  is  a  function  of  potential  on  charge,  i.e. 
this  ratio  increases  towards  completion  of  charging  which  is 
not  normally  the  case  for  ideal  battery  cells. 

A  maximum  in  PD  with  respect  to  ED  at  a  particular  SOC 
defined  by  an  open-circuit  potential  is  also  predicted  the¬ 
oretically  for  capacitors.  Given  that  for  a  particular  SOC  of  a 
capacitor 

V=V,-/R 
on  discharge 
P=i{V)=iV-i2R 
and 

E=jCV2=^iV-iR)2 

the  slope  of  the  ln-ln  Ragone  relation  is,  therefore 
dlnP  E dP/di _l/2CV2(V-2iR)  Q 
d  In  E  PdE/di  iV(-CRV) 


The  ED  corresponding  to  the  maximum  PD,  viz.  = 
V,  2/4 R  [10],  corresponding  to  the  conditions  given  by  Eq. 
(29),  is  then 


3.  Results  and  discussion 

i.  /.  Quantitative  evaluations  of  polarization  effects  in 
battery  discharge 

In  Section  2.2  we  have  given  equations  which  represent 
the  effects  of  (i)  ohmic  (iR)  polarization  and  (ii)  kinetic 
(Tafel)  polarization  on  both  the  ED  and  the  PD  of  battery 
systems  on  discharge  which  characterize  the  kagone  plot  for 
a  given  system.  Here  we  show  by  means  of  some  quantitative 
numerical  calculations  how  the  following  parameters  of  bat¬ 
tery  systems  determine  the  shape  and  position  of  Ragone  plots 
between  the  scales  of  ED  and  PD.  Both  these  quantities 
depend  (cf.  Eqs.  (6)  and  ( 18))  on  the  electric  polarizations 
suffered  at  the  battery  electrodes  which  increase  logarithmi¬ 
cally  with  i  in  the  kinetic  (Tafel)  regime  and  linearly  in  i  in 
the  ohmic  (iR)  regime.  Normally  both  factors  have  to  be 
considered  except  at  extremes  of  discharge-current  values. 

The  parameters  are  ( i)  the  internal  ohmic  resistance,  R,  of 
the  battery  system  (electrolyte  and  electrode  structures)  and 
(ii)  the  Tafel  equation  parameters  a  (or  r0)  and  b  character¬ 
izing  the  dependence  of  i)a  on  In  /.  Both  the  R  and  the  Tafel 
parameters  have  been  recognized  qualitatively  as  determining 
factors  in  Ragone  plots  [  1J,  e.g.  for  fuel-cell  performance 
[  1 1  ] ,  but  quantitative  numerical  evaluations  for  battery  sys¬ 
tems  have  been  lacking,  especially  the  separate  roles  of  the 
values  of  the  Tafel  slopes  for  the  anode  and  cathode  reactions. 
Ragone  plots  themselves  are  usually  shown  on  an  empirical 
basis  with  little  or  no  reference  to  the  relative  importance  of 
the  ohmic  and  kinetic  contributions  to  the  PD  versus  ED 
relations.  It  is  the  aim  of  the  ‘theoretical'  calculations,  the 
results  of  which  are  shown  diagrammatically  in  the  figures 
below,  to  demonstrate  the  relative  influences  of  the  ohmic 
and  kinetic  effects  for  arbitrary  but  realistic  selections  of 
parameter  values.  Although  the  parameters  we  have  chosen 
to  represent  Ragone  plots  for  various  conditions  are  arbitrary, 
they  have  been  selected  to  correspond  to  ranges  of  practical 
behaviour  commonly  observed.  In  the  case  of  the  ohmic 
polarization  factor,  which  diminishes  cell  voltage  and  hence 
both  PD  and  ED,  it  is  to  be  noted  that  i  x/f  scale  together  so 
that  the  equivalent  behaviour  arises,  for  example,  with  high 
i  and  low  R  as  for  low  i  and  high  R.  The  important  factor 
deteimining  the  shape  and  position  of  Ragone  plots  in  the  PD 
versus  ED  plane  is  how  the  iR  polarization  compares  with 
increasing  Tafel  polarization  as  i  is  increased  with  increasing 
operating  power,  see  Eqs.  (6)  or  (7).  Also,  the  separate 
dependences  of  the  ED  and  PD  on  i  are  useful  to  evaluate  in 
relation  to  corresponding  Ragone  plots. 
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The  following  are  the  bases  of  the  calculations  for  simu¬ 
lated  Ragone  plots  for  batteries: 

(i)  An  overall  polarization  equation  for  electrode  potential 
Vas  a  function  of  i  (including  iR  and  17,,,  but  excluding  ifc.) 
is  written  (cf.  Eq.  (6))  as 

V=Al'rev-;/?-(a  +  £logz)  (30) 

Nominal  values  of  AVrcv=  1  V  and  i0=  10“ 5  A  kg"1  arc 
taken.  The  Tafel  constant  a  is  —  b  log  /0. 

(ii)  Energy-density,  E,  is  taken  as  <2  X  V  where  Q  is  the 
charge  capacity  in  MC  kg~ :  (about  10  F),  nominally  taken 
as  1  for  the  purpose  of  calculations  of  relative  behaviours. 

(iii)  Power-density,  Py  is  taken  as  i  X  V. 

( iv)  V,E  and  P  were  calculated  for  R  =  0.01 , 0. 1 , 1 .0  and 
10  12,  and  values  of  the  Tafel  parameters  a,  b,  respectively, 
as  0.6  and  0.120  V,  0.21  and  0.042  V,  0.15  and  0.03  V  and 
limitingly  as  0  and  0  V  ( no  kinetic  polarization) .  Other  values 
of  a,  corresponding  to  various  other  selected  values  of  In  i0, 
could,  of  course,  be  taken. 

The  calculations  resulted  in  a  family  of  16  log  P  versus  log 
E  curves;  the  range  of  i  values  for  each  curve  was  selected 
such  that  V>0,  i.e.  P>  0,  £>0,  and  that  the  curves  covered 
the  range  of  i  values  within  which  the  expected  maximum  in 
P  (Eq.  (10))  would  be  demonstrated.  Note  that  Eq.  (30) 
predicts  (when  multiplied  by  i )  an  initial  increase  in  P  with 
increasing  i  when  i  is  sufficiently  small  that  the  two  polari¬ 


zation  terms  which  diminish  V  from  A  Vrov  are  negligible  but 
a  decrease  in  P  at  larger  i  values;  hence  a  maximum  in  P 
arises,  as  was  shown  analytically  in  Eqs.  (9),  (10)  and  (12). 
However,  it  is  interesting  that  almost  all  empirically  deter¬ 
mined  Ragone  plots  (e.g.  Fig.  I )  do  not  exhibit  maxima  but 
only  a  roll-off  of  P  with  increasing  i  and  hence  diminishing 
E.  This  is  presumably  because,  with  practical  systems,  suf¬ 
ficiently  high  discharge  currents  were  not  used  in  the  per¬ 
formance  tests  to  reach  and  exceed  the  maximum  in  P,  or 
alternatively  current  limitation  by  diffusion  control  had  set 
in,  as  may  arise  in  fuel  cells  and  some  batteries,  e.g.  due  to 
consumption  of  Li +  ions  or  to  pH  changes. 

The  following  summarizes  the  set  of  log  P  versus  log  £, 
and  P  or  E  versus  log  i  relations  calculated: 

(i)  Fig.  2(a)-(d),  constant  b  value,  various  values  of 
R  =  0.01, 0.1, 1.0,  10.012; 

(ii)  Fig.  3(a)-(d),  constant  R,  various  values  of 
b =0.120, 0.042,  and  0.030  V; 

(iii)  Fig.  4,  constant  b =0. 120  V,  various  values  of  R  for 
P  and  for  £  as  a  function  of  log  i;  and 

(iv)  Fig.  5,  R=  1  n,  with  various  b  values  for  log  P  and 
for  E  versus  log  i. 

It  wii!  ^  useful  now  to  present  a  commentary,  in  a  com¬ 
parative  way,  on  the  forms  of  these  results  as  shown  in  the 
above  figures.  In  all  the  computed  series  of  Ragone  plots  in 
Figs.  2(a)-(d)  and  3(a)-(d),  ED  at  first  decreases  almost 


Fig  2.  Calculated  Ragone  diagrams  for  a  battery  having  A  V(cv  =  l  V,  i„=  10"5  A  kg-  \  showing  the  effect  of  ohmic  polarization  for  fixed  Tafel  polarization: 
(a)  b= 0. 120  V.  a  -  0.60  V;(b)f»  =  0.042  V.  a  =  0.21  V;  (c)  6  =  0.030  V.c  =0.15  V;  <d)  fc~0V.a  =  0V. 
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En&ffly-density  x  Q>  (V)  Energy-density  x  CP  (V) 

Fig,  3.  Calculated  Ragone  diagrams  for  a  battery  having  A  Vrev=  1  V.  i'0=  10”  5  A  kg- showing  the  effect  of  Tafel  polarization  for  fixed  ohmic  polarization: 
(a)  R =0.01  ft;  (b)  /f =0.10  ft;  (c)  /?=  1.0ft,  and  (d)  /f=  10.0  ft. 


Current-density  (A  kg*’) 

Fig.  4.  The  effect  of  ohmic  polarization  on  PD  and  reduced  ED  as  a  function  of  i  for  a  battery  characterized  bj  A  Vm=  1  V.  <<,=  10-5  A  kg-1, 6=0.120  V, 
and  o=0.60  V. 

linearly  with  increasing  PD  (on  the  log-log  plots)  in  the  first  (as  expected  according  to  Eq.  ( 12)  ).  The  roll-off,  before  the 

decade  of  values  of  ED  0.1  to  1  Wh  kg  ”1.  Then  the  PD  maximum  power  is  attained,  takes  place  at  lower  power  as  R 

maximum  (Eqs.  (9),  (10)  and  (12))  is  attained  followed  is  increased  from  0.01  to  10  ft  or  more  as,  is  intuitively 

by  progressively  decreasing  PD  with  decreasing  ED,  the  PD  expected. 

levels  attained  being  progressively  lower  (as  expected)  for  Fig.  2(a)  is  calculated  for  b= 0. 120  V,  a = 0-6,  Fig.  2(b) 

the  larger  R  values  (greater  iff-drop  component  values).  The  for  b =0.042  V,a= 0.21,  Fig.  2(c)  forh  =  0.029  V,a=0.15 

slopes  of  the  log  P  versus  log  E  lines  below  the  Pmtx  are  and  Fig.  2(d)  (reference  case),  korfr=0,n  =  0  (zero  kinetic 

approximately  the  same  for  small  i  and,  as  i ->0  approach  1  polarization).  Fig.  2(d)  thus  illustrates  the  limiting  case 
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Fig.  5.  The  effect  of  Tafel  polarization  on  PD  and  reduced  ED  as  a  function  of  i  for  a  battery  characterized  by  A  Vrev  =  1  V,  io=  10"5  A  kg  -  l ,  and  R  =  I  ft. 


where  only  i7?  polarization  determines  the  shape  of  the 
Ragone  plots.  The  power,  is  of  course,  attenuated  as  R  is 
increased  from  0.01  to  10  ii,  or  more. 

The  b  values  were  chosen  as  ;hoKr  typically  considered 
and  often  observed  for  the  hydrogen-evolution  reaction  and 
other  well-known  electrode  processes  {7,12]. 

The  series  of  Fig.  3(a)-(d)  shows  Ragone  plots  resulting 
from  variation  of  Tafel  slope  b  values  with  a  selected  fixed  R 
value  taking  R  =  0.01 ,0.1,  1.0  and  10  ft  each  for  a  series  of 
three  b  values. 

The  same  general  forms  of  the  resulting  Ragone  plots  in 
the  series  ‘3’  calculations  arise  as  in  the  series  ‘2’  but  the 
variation  of  b  mainly  affects  the  shape  of  the  Ragone  plot  at 
the  high  ED  (low  i)  end  before  the  maximum  in  P  arises 
( with  decreasing  E) .  The  curves  again  have  the  same  slopes 
(parallel  lines)  beyond  the  maxima,  at  decreasing  E  values, 
almost  independently  of  the  assigned  Tafel  slope  value. 
Again,  the  power  levels  attained  diminish  with  increasing 
values  of  the  chosen  R,  0.01  to  10  ft.  It  is  the  high  ED  ends 
of  the  Ragone  plots  that  are  most  affected  by  the  assigned  b 
value,  i.e.  where  the  kinetic  polarization  increases  with 
increasing  f  the  more  rapidly  the  greater  is  the  b  value  in  the 
series  b = 0.029  <  0.042  <  0. 1 20  V. 

3.2.  Individual  dependences  of  P  and  E  on  log  i 

In  the  evaluation  of  battery  power  systems,  the  Ragone 
relations  are  not  usually  expressed  complementarily  in  terms 
of  the  individual  P  and  E  values  as  a  function  of  i.  However, 
this  is  an  informative  way  of  examining  the  data  for  a  given 
power  source.  Using  the  quantities  calculated  for  Figs.  2  and 
3,  Figs.  4  and  5  show  the  dependences  of  both  P  and  the 
reduced  ED  (E/Q)  on  log  i  for  variation  of  R  from  0.01  to 
10 ft,  for  b  =  0. 120  V  ( Vy  =  1 .0  V)  (Fig.  5)  and  in  Fig.  6  for 
variation  of  b  taken  as  0.029,  0.042,  0.120  V  with  a  single 
value  for  R,  taken  as  1  ft.  These  show  clearly  how  the  PD 
falls  off,  after  the  maximum,  only  strongly  at  high  i  in  the 
decade  0.1  to  1  A  kg' There  is  a  broader  fall-off  of  E/Q 


both  with  increasing  R  (0.01  to  10ft,  Fig.  5)  and  increasing 
b  (Fig.  6).  The  plots  for  ED  versus  log  i  for  increasing  R 
closely  resemble  the  form  of  Tafel  relations  for  a  given  a  ( or 
log  i0)  and  b  value  with  increasing  solution  resistance  in  series 
with  the  electrode. 

Finally,  in  this  analysis  of  the  relation  between  effects  of 
iR  and  7ja  polarization,  it  is  informative  to  show  the  progres¬ 
sive  changes  of  the  form  of  a  ‘Tafel’  plot  with  increasing  iR- 
drop  component  in  the  overall  iR  +  t/.,  polarization.  This  can 
be  done  in  a  generalized  way  by  plotting  reduced  current 
density  as  i/i0  on  the  axis  of  abscissae  agai  nst  a  polarization 
quantity  rj  =  i?a  +  iR  with  iR  taken  as  i()R  X  i/i0  for  a  series  of 
values  of  R  expressed  as  /„/? = 0  to  IqR  —  2  X  1 0  “  2,  as  shown 
in  Fig.  6.  Here  the  transition  from  pure  Tafel  (logarithmic) 
polarization  behaviour  (curved  on  the  linear  plot  versus  i/i0 
with  i||ff-0)  to  much  more  linear  behaviour  (on  the  linear 


Fig.  6.  Effect  of  combination  of  a  linear  (ohmic)  and  exponential  (Tafel) 
polarization  effect  on  generalized  cunent  vs.  potential  curves  for  an  electrode 
process.  Current-density  scale:  relative  i/i„;  overpotential  scale:  17  + iff 
(iff=/yff  Xi/i0).  Diagram  shows  transition  from  mainly  ‘Tafel’  to  mainly 
‘ohmic’  polarization  for  a  Tafel  slope  of  0. 1 18  V  and  i = i<,x  I0',,b. 
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i//0plot)  is  seen  with  increasing  values  ofthcroR  polarization 
component. 

3.3.  Ohmic  and  Tafei  polarization  components  of  Ragone 
plots 

it  is  interesting  and  useful,  for  comparative  practical  eval¬ 
uation  of  electrochemical  power  sources,  to  evaluate  and  thus 
distinguish  separately  the  ohmic  and  the  Tafei  polarization 
components  that  determine  the  form  of  Ragone  plots. 

In  Fig.  7  are  plotted  reduced  EDs,  i.e.  E/Q ,  a  quantity  in 
V,  against  log  r.  Fig.  7  shows  the  components  of  this  voltage 
fo  a  given  Tafei  slope  {b  =  0. 120  V)  and  a  variable  iR com¬ 
ponent  with  R  varying  from  0.0!  to  10  ft,  line  A  shows  the 
reduced  ED  prior  to  polarization  corrections,  line  B  represents 
the  ED  losses  due  to  Tafei  polarization,  lines  Cl ,  C2,  C3  and 
C4  the  ED  losses  due  to  ohmic  effects,  and  lines  El,  E2,  E3 
and  E4  the  net  ED  of  the  systems.  This  figure  can  be  analysed 
:o  determine  the  i  range  over  which  each  of  the  polarization 
terms  dominates  the  behaviour  of  the  resultant  ED.  At  low  i, 
the  ED  is  linear  in  log  i,  and  clearly,  the  Tafei  polarization 
component,  which  is  much  larger  than  the  iR  component  in 
this  region,  is  dominant.  For  a  fixed  Tafei  slope,  and  variable 
R,  the  current  range  over  which  the  ED  versus  log  i  curve  is 
linear,  increases  as  R  decreases.  At  high  current  densities,  the 
iR  polarization  increases  rapidly  (cf.  Fig.  7)  relative  to  the 
Tafei  polarization  component  and  results  in  a  sharp  decline 
in  the  ED.  A  similar  analysis  of  E/Q  for  fixed  R  and  varied 
b  shows  that  a  larger  value  of  b  leads  to  a  greater  component 
of  kinetic  polarization  for  a  given  value  of  other  things 
being  equal.  Thus,  the  diminution  of  ED  can  be  quite  signif¬ 
icant  for  large  b  even  at  small  i.  This  component  analysis 
shows  that  r)u  is  most  important  at  high  ED  (low  /)  rather 
than  at  high  PD  (higher  i). 

Fig.  8  shows  plots  of  the  contributions  of  AVrev,  iR  and 
(a  +  b  log  i)  to  the  PD.  Each  of  the  terms  in  Eq.  ( 6) ,  rA  V^, 


i2R  and  i(a  +  b  log  i)  were  calculated  for  specific  parameter 
values  (/?,  a,  and  b)  and  are  plotted  versus  log  i.  The  curves 
in  Fig.  8  are  for  a  fixed  iR  component  {R=  1.0  ft)  and  a 
varied  Tafei  slope,  b  —  0. 1 20, 0.029, 0  V  and  have  been  plot¬ 
ted  as  a  linear  log  relation,  as  for  the  E/Q  diagram,  line  A  is 
the  polarization-free  PD,  B  the  iR  polarization  contribution. 
Cl  and  C2,  the  Tja  components;  PI,  P2,  and  P3  are  the  polar¬ 
ization-corrected  PD’s. 

The  interpretation  of  these  figures  is  more  difficult  than  for 
the  complementary  ED  diagrams  previously  discussed  for 
which  the  reduced  ED,  before  correction  for  polarization 
losses,  was  constant  for  all  i's  (Fig.  7,  line  A  =  1 ).  The  PD 
before  correction  for  polarization  losses,  is,  of  course,  directly 
dependent  on  i  as  defined  by  i  X  A  Vrev.  Furthermore,  neither 
the  decrease  in  PD  due  to  7}a  nor  to  iR  is  linear  in  log  i. 

Several  general  statements  regarding  these  diagrams  can 
be  made;  again,  as  for  the  case  of  ED,  the  activation  polari¬ 
zation  is  of  greater  significance  at  low  i,  and  the  iR  contri¬ 
bution  is  larger  at  high  i.  Interestingly,  it  is  only  when  both 
irja  and  i-R  are  significant  with  respect  to  iA  V„.v  that  maxima 
in  the  PD  curves  arise. 

At  sufficiently  low  /,  ED  is  significantly  reduced  as  i 
increases  (Fig.  7)  and  is  controlled  by  rja,  whereas  PD  is 
relatively  independent  of  i  (Fig.  8)  and  is  relatively  small. 
At  intermediate  i,  ED  decreases  more  rapidly  as  both  iR  and 
i?a  reduce  the  operating  voltage,  (Fig.  7).  Power  density 
increases  as  the  iA  Vrcv  term  dominates  the  i2R  and  i%  polar¬ 
izations.  At  high  i,  ED  drops  sharply  towards  zero  and  PD 
reaches  a  maximum  and  then  decreases.  Energy-density 
decreases  over  the  full  range  of  /,  whereas  the  change  in  PD 
is  significant  only  at  intermediate  and  high  i  values. 

3.4.  Discussion  of  calculations  for  power-density  versus 
energy-density  of  a  capacitor 

The  following  are  the  bases  of  the  calculation  s  for  simu¬ 
lation  of  Ragone  plots  for  capacitors: 


Fig.  7.  Contribution  of  ohmic  and  Tafei  polarization  to  reduced  ED  as  a  function  of  i  fora  battery  characterized  by  A  V^.v  =  !  V,  i0=  IQ-5  A  kg-1.  6  =  0.120 
V,  and  a  =  0.60  V.  Line  A  corresponds  to  the  ED  of  a  battery  prior  to  polarization  corrections,  line  B  represents  ED  losses,  due  to  Tafei  polarization,  lines  Cl. 
C2,  C3  and  C4  to  ED  losses  due  to  ohmic  effects  for  R  =  10.  1.0. 0.1,  and  0.01  £1.  respectively  and  lines  El,  E2,  E3  and  E4  the  net  ED  of  the  systems. 
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Current-density  (A  kg  ') 

Fig.  8.  Contributions  of  ohmic  andTafe!  polarization  to  PD  as  a  function  of  i  for  a  battery  characterized  by  A  Vrey  =  I  V,  /„=  10" s  A  kg'1,  /?=  1.0  ft.  Line  A 
corresponds  to  the  PD  of  a  battery  prior  to  polarization  corrections.  line  B  represents  PD  losses  due  to  ohmic  polarization,  lines  C  t  and  C2  to  PD  losses  due  to 
Tafel  effects  for  fc=0.  !20  and  0.030  V,  respectively,  lines  PI.  P2  and  P3  the  net  PD  of  the  systems. 


( i)  An  overall  polarization  equation  for  electrode  potential 
as  a  function  of  i  is  taken  as  V- V, ■-  iR. 

b)  Energy  -  densi  ty,  £,  is  taken  as  \/2QV=  M2CV2,  where 
Q  is  the  charge  capacity  in  MC  kg" ',  nominally  taken  as  1 
MC  kg"  1  (about  10 F). 

(iii)  Power-density,  P,  is  taken  as  /X  V. 

Fig.  9  illustrates  the  Ragone  plots  for  fully  charged  capac¬ 
itors  having  constant  Q  (  —  1  MC  kg" 1 ),  and  variable  nom¬ 
inal  C  (  =2.0,  J.0,  0.5  and  0.1  MF  kg"1)  values.  Each  of 
these  cases  defines  an  initial  open-circuit  potential  for  a 
charged  capacitor  of  (  =0.5,  1,  2,  and  10  V).  The  effect 
of  iR  polarization  is  included  in  this  analysis  for  R  =  0. 1 ,  1.0 
and  10  mil.  As  discussed  by  Miller  [10],  a  capacitor  is 
considered  effectively  discharged  when  the  operating  poten¬ 
tial  falls  to  VJ2 ,  at  which  time  the  capacitor  has  delivered 
75%  of  its  charge  capacity  into  a  load  resistance.  For  a  con¬ 
stant-load  discharge,  the  remaining  25%  of  capacity  can  only 
be  recovered  over  an  extended  time  and  increasingly  lower  i 


(remembering  that  current  decreases  exponentially  in  time 
through  a  constant  resistance  load).  Thus  the  curves  in  Fig. 
9  correspond  to  potentials  between  V,  and  VJ2.  Decreasing 
C  while  maintaining  Q  constant,  increases  the  initial  operat¬ 
ing  voltage  and  thus  both  the  maximum  ED  and  PD  attainable. 
If  the  capacitor  is  discharged  past  V,/2,  both  ED  and  PD  will 
decrease,  cf.  Eq.  (29).  Obviously,  decreasing  the  iR  polari¬ 
zation  increases  the  maximum  PD.  The  maximum  ED  is 
effectively  independent  of  R ,  as  it  always  arises  when  i  0. 

Fig.  10  shows  how  the  SOC  dependence  of  ED  and  PD 
values  can  be  calculated  as  a  function  of  i.  Here  0,  (Eq. 
(22))  is  taken  as  1  MC  kg" l,  Cas  1  MF  kg"  \  as  1  V  and 
R  as  10  mil.  The  resulting  ED/PD  relationships  were  cal¬ 
culated  as  a  function  of  i  for  a  capacitor  holding  56,  36  and 
25%  of  the  charge  capacity  (on  an  Ah  basis)  of  the  fully 
charged  capacitor.  Clearly  (sec  Fig.  10)  both  ED  and  PD  are 
dependent  on  the  SOC  of  the  capacitor,  and  the  maximum 
PD  deliverable  by  the  capacitor  decreases  as  the  capacitor  is 


Energy-density  (kW  ks  kg') 

Fig.  9.  Calculated  Ragone  plots  for  fully  charged  capacitors  having  constant  Q=  1  MC  kg“  \  variable  nominal  C  ( =2.0.  1.0, 0.5  and  0.1  MF  kg1)  values, 
and  variable  ohinic  polarization  (R  =  0. 1 .  1.0,  1 0.0  mil } 
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Energy-density  (kW  ks  kg*1) 

Fig.  10.  Calculated  Ragone  plot  showing  the  SOC  dependence  of  ED  and  PD  for  a  capacitor  characterized  by  Q,~  1  MC  kg"  \  C  as  I  MF  kg- 1.  V,  =  1  V  and 
/?=*  10.0  mil  holding  100,  56, 36  and  25%  of  the  charge  capacity  (on  Ah  basis)  of  the  fully  charged  capacitor. 


discharged.  This  is  a  significant  practical  characteristic  that 
must  be  considered  during  system  design  of  battery /capacitor 
hybrids. 

Fig.  1 1  illustrates  normalized  ED  and  PD  versus  i  plotted 
as  a  function  of  the  SOC  for  conditions  as  described  for  Fig. 
10.  Fig.  1 1  reinforces  the  significance  of  the  dependence  of 
both  ED  and  PD  on  the  SOC  of  the  capacitor,  and  illustrates 
that  there  is  a  maximum  power  point  (MPP)  of  operation 
and  that  the  MPP  is  dependent  on  SOC.  The  MPP  corresponds 
to  the  current  density  and  device  potential  at  which  maximum 
power-density  may  be  delivered.  In  Fig.  II,  the  dark  solid 
lines  correspond  to  the  EDs  and  PDs  up  to  the  cutoff  voltage 
for  a  discharged  capacitor  of  Vj2\  the  dashed  lines  to  an  ED 
limit  of  0.  As  the  SOC  is  decreased,  the  energy  limit  (viz. 
25%  of  the  energy  of  the  fully  charged  capacitor  remaining 


at  the  cutoff  voltage  of  V,/2)  is  reached  at  progressively  lower 
i’s,  and  clearly  this  limit  is  attained  at  /* s  smaller  than  those 
defined  by  the  MPP.  It  is,  therefore,  impossible  to  operate  at 
the  MPP  for  any  SOC  other  than  100%,  illustrating  that  the 
choice  of  the  optimum  capacitor  (i.e.  selection  of  £),  C,  V#) 
for  a  particular  application  is  not  a  simple  problem. 

3.5.  Kinetic  polarization  factor  in  Ragone  plots  for  a 
pseudocapacitance  device 

In  Section  3.4  the  polarization  of  the  capacitor  was 
assumed  to  be  ohmic  ( V=  Vt  -  iR)  due  to  the  internal  resis¬ 
tance  of  the  porous  electrode  matrix  and  of  the  electrolyte 
within  it.  Such  conditions  apply  principally  to  the  ‘double- 
layer’  type  of  supercapacitor.  However,  when  the  charge  stor- 


Current-density  (A  kg-*) 

Fig.  11.  Normalized  ED  and  PD  vs.  i  ptoued  as  a  function  of  SOC  for  a  capacitor  characterized  by  Q,=  1  MC  kg-  \  Cas  l  MF  kg-*,  I  V  and  /?=  10.0 
ran  holding  100,  56  and  36%  of  the  charge  capacity  (on  Ah  basis)  of  the  fully  charged  capacitor.  Solid  lines  correspond  to  the  EDs  and  PDs  up  to  the  cutoff 
voltage  for  a  discharged  capacitor  of  V,/2;  the  dashed  lines  to  an  ED  limit  of  0. 
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age  of  the  device  originates  principally  from  faradaic 
pseudo-capacitance  (e.g.  associated  with  electrosorption 
( 13]  or  solid-state  or  surface  redox  processes,  as  with  Ru02 
19]),  then  the  analysis  for  battery-type  polarization  relations 
developed  in  Section  2.2  will  apply,  i.e.  when  both  the  ohmic 
and  Tafel  polarization  factors  are  involved. 

However,  in  the  case  of  a  supercapacitor  energy-storage 
device,  independently  of  whether  the  capacitance  arises  elec¬ 
trostatically  from  the  interfacial  double-layer  or  faradaicaliy 
mainly  from  pseudo-capacitance,  the  potential  will  decline 
with  extent  of  discharge  uniformly  if  the  capacitance  were 
constant  or  in  a  more  complex  way  if  the  capacitance  itself 
has  some  dependence  on  potential  as  is  often  the  case  exper¬ 
imentally  (9,13,14]. 

This  factor  of  decline  in  potential  with  extent  of  discharge 
is  a  general  and  intrinsic  characteristic  of  any  type  of  capacitor 
device;  it  is  the  main  factor  which  distinguishes  the  behaviour 
of  batteries  from  that  of  capacitors  reflected  in  the  qualitative 
difference  between  Ragone  plots  for  capacitors  and  those  for 
batteries,  with  the  exception  (cf.  Ref.  [8])  of  the  lithium- 
intercalation  type.  This  difference  is  illustrated  in  the  present 
paper  through  Figs.  9-1 1  in  relation  to  Ragone  plots  for 
battery-type  energy  storage  systems  shown  earlier. 


4.  Conclusions:  interpretation  of  practical  Ragone 

diagrams 

In  Figs.  3  and  4  we  showed  that  the  maximum  in  PD  for  a 
battery  system  defined  by  an  open-circuit  potential  of  V%  is 
primarily  determined  by  the  iR  losses,  provided  that  R  is  not 
negligibly  small,  and  that  the  negative  slope  of  PD  versus  ED 
is  determined  by  losses  associated  with  the  tju.  A  large  neg¬ 
ative  slope  indicates  a  low  rja,  and  a  smaller  slope,  a  high 
Obviously,  a  low  and  therefore  a  sharp  decrease  in  PD  at 
high  ED  is  preferred.  The  maximum  in  ED  arises  as  /-*0, 
where  losses  associated  with  i?a  are  significantly  ,r>rger  than 
iR  losses.  Thus,  a  low  r)a  also  favours  a  high  ED.  \  ' 

From  Figs.  8-10  it  is  seen  that  the  Ragone  plot  for  a  fully 
charged  capacitor  is  effectively  at  a  constant  PD  over  most 
of  the  range  of  ED.  Further,  and  of  major  practical  impor¬ 
tance,  such  a  Ragone  plot  is  substantially  dependent  on  the 
SOC,  and.  as  the  capacitor  discharges,  the  log  P  versus  log  E 
dependence  becomes  more  sloping  (Fig.  9).  Matching  of 
capacitor/batlery  systems  should  take  into  account  this 
dependence  of  the  capacitor  PD  versus  ED  relation  on  SOC. 


Thus,  in  an  hybrid  capacitor  load-levelled  battery  system  for 
electric-vehicle  application  [3,4],  optimal  performance  of 
ciie  battery-capacitor  combination  will  be  achieved  along  the 
initial,  almost  flat  part  of  the  Ragone  curves  where  PD 
remains  almost  constant.  Obviously  also,  as  seen  from  Fig. 
10,  the  best  power  performance  provided  by  the  capacitor 
component  will  be  realized  when  the  system  is  operated  at 
the  highest  possible  SOC,  e.g.  maintained  by  recharge  from 
the  battery  component  on  level  and  uniform-speed  driving 
and/or  from  a  regenerative  braking  system. 
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